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Abstract 

The cat-eyes steady state solution in the framework of hydrodynamics 
describing an infinite row of identical vortices is extended to the magnetohy- 
drodynamic (MHD) equilibrium equation with incompressible flow of arbi- 
trary direction. The extended solution covers a variety of equilibria including 
one- and two-dimensional generalized force-free and Harris-sheet configura- 
tions which are preferable from those usually employed as initial states in 
reconnection studies. Although the vortex shape is not affected by the mag- 
netic field, the flow in conjunction with the equilibrium nonlinearity has a 
strong impact on isobaric surfaces by forming pressure islands located within 
the cat-eyes vortices. More importantly, a magnetic-field-aligned flow of ex- 
perimental fusion relevance and the flow shear have significant stabilizing 
effects in the region of pressure islands. The stable region is enhanced by an 
external axial ("toroidal") magnetic field. 
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I. Introduction 



Sheared flows influence the equilibrium and stability properties of mag- 
netically confined plasmas and result in transitions to improved modes either 
in the edge region (low-to-high-mode transition) or in the central region (in- 
ternal transport barriers) of fusion devices. As concerns equilibrium, the 
convective velocity term in the momentum equation makes the isobaric sur- 
faces to deviate from magnetic surfaces, unlike the case of quasistatic steady 
states [1], thus potentially affecting stability. For symmetric two dimensional 
equilibria, this effect has been examined on the basis of analytic solutions to 
linearized forms of generalized Grad-Shafranov equations, e.g. Ref. [2]. For 
flows of fusion concern, i.e. for Alfven Mach numbers of the order of 0.01, 
this deviation is small and consequently isobaric and magnetic surfaces have 
the same topology. 

Aim of the present study is to examine the impact of flow in conjunction 
with nonlinearity to certain equilibrium and stability properties in relation 
to the departure of the isobaric from magnetic surfaces. Motivation was a 
solution of a nonlinear form of the hydrodynamic equation describing the 
steady motion of an inviscid incompressible fluid in two dimensional plane 
geometry, known as "cat eyes" , which represents an infinite row of identical 
vortices ([3111]; see also Fig. [I]). This solution is extended here to the MHD 
equilibrium equation with incompressible flow. Then, the stability of the 
extended solution is examined by means of a recent sufficient condition [5]. 
The major conclusion is that owing to the nonlinearity of the equilibrium, the 
flow and flow shear affect drastically the pressure surfaces and have significant 
stabilizing effects in the region of modified pressure. 

The MHD equilibrium equations with incompressible flow for transla- 
tionally symmetric plasmas are reviewed in Sec. II. In Sec. Ill a solution 
of the pertinent generalized Grad-Shafranov equation describing a whole set 
of equilibria is constructed as an extension of the cat-eyes solution. Then 
for parallel flows and constant density the stability of the solution obtained 
is studied in Sec. IV. Sec. V recapitulates the study and summarizes the 
conclusions. 

II. Review of the equilibrium equations 

The MHD equilibrium states of a translationally symmetric magnetized 
plasma with incompressible flows satisfy the generalized Grad-Shafranov 
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equation [5|, 



;i - M 2 )v 2 >ip - -(M-yiv^r + Ps + 



p? 



(i) 



and the Bernoulli relation for the pressure 

p = p s ^)-\m 2 \v^ 



(2) 



Here, the function ip(x, y) labels the magnetic surfaces with (x, y, z) Cartesian 
coordinates so that z corresponds to the axis of symmetry and (x, y) are 
associated with the poloidal plane; M(ip) is the Mach function of the poloidal 
velocity with respect to the poloidal-magnetic-field Alfven velocity; B z is the 
axial magnetic field; for vanishing flow the surface function P s (ip) coincides 
with the pressure; the prime denotes a derivative with respect to ip. The 
surface quantities M{jp), B z {jp) and P s (ip) are free functions for each choice 
of which ([1]) is fully determined and can be solved whence the boundary 
condition for ip is given. Also, to completely determine the equilibrium, a 
couple of additional surface functions are needed, i.e, the density, q(i/j), and 
the electrostatic potential, $(V0- Details including derivation of (Op) and ([2]) 
can be found in Refs. [61 [7] . 

Eq. (CQ) can be simplified by the transformation [HI [9] 



u(ip) 



r l - M 2 (g)] dg 



which reduces flU) to 



0. 



(3) 



(4) 



Also, (j2J) is put in the form 

P = PJu) 



M 2 



2(1 — M 2 ) 



|Vn| 



(5) 



Note that (j3j) free of a quadratic term as |Vu| 2 is identical in form with the 
quasistatic MHD equilibrium equation as well as to the equation governing 
the steady motion of an inviscid incompressible fluid in the framework of 
hydrodynamics. Transformation Q does not affect the magnetic surfaces, 
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it just relabels them. Also, once a solution of fll]) is found, the equilibrium 
can be completely constructed in the w-space; in particular, the magnetic 
field, current density, velocity, and electric field can be determined by the 
relations: 



B 
j 



B z e z + ^ 

1 — M 2 
dB 



- M 2 

-1/2 



-1/2 



V z u + -- 



e z xVu 

ldM 2 



2 du 



M' 



-3/2 



E 



du 
— B 

— -Vu. 

du 



■e 7 xVu 



1/2 d$ 

du 



(6) 
(7) 

(8) 
(9) 



Analytic solutions to linearized forms of (J4]) have been constructed for 
quasistatic [TDl E] and stationary equilibria [TJ. As already mentioned in 
Sec. I, for flows of experimental fusion relevance (|M| ~ 0.01) the departure 
of the isobaric from magnetic surfaces is small (see for example Fig. 2 of Ref. 
[2]), so that the topology of these two families of surfaces is identical. 



II. Magnetohydrodynamic "cat eyes" with flow 



The present section aims at extending the hydrodynamic cat-eyes solution 
to (jl]) and examine certain equilibrium characteristics in connection with the 
impact of the flow together with nonlinearity. For convenience we introduce 
dimensional quantities: x = x/L, y — y/L, u — u/(B z0 L), g = g/ go, P = 
P/B 2 z0 , B = B/B z0 , j = j/(B z0 /L), v = v/v A0 , where v A0 = B^/^/g^ , and 
E = E/(i?^ f J 4 ); here, L, B z0 , and £ are reference quantities to be defined 
later. Eqs. (j!J) and (jSJ) hold in identical forms for the tilted quantities and 
will be further employed as dimensionless by dropping for simplicity the tilde. 
To construct a cat-eyes solution we make the ansatz 

^H.'-D^-H do) 

du 

by which (JIJ) reduces to the following form of Liouville' s equation: 

V 2 n = (1 - e 2 )exp(-2n). (11) 
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Eq. (TTTT) admits the solution 



w = In [cosh(y) — ecos(x)] , (12) 

the characteristic lines of which are shown in Fig. [TJ The parameter e 
determines the vortex size; for e = 1 the solution represents an infinite row 
of point vortices and for e = it becomes one-dimensional: u = In cosh y. 
It is noted here that though ffl2|) is singular in the limit of y — > oo, all the 
local equilibrium quantities are everywhere regular. Eq. ([TO]) can be solved 
for P s (u) + B 2 z /2 to yield 

B 2 1 — e 2 1 - e 2 

Ps + -f = exp (-2u) + c = — - — -2 + c , (13) 

* * 2 (cosh y — e cos 

where c is a constant. The equilibria described by (fT2l) and (|T3|) have the 
following characteristics: 

1. The vortices are by construction of solution ( Tl2l) identical to the re- 
spective hydrodynamic vortices, viz., the magnetic field does not affect 
the vortex shape. 

2. Since magnetic field and current density lie on the velocity or magnetic 
surfaces, the vortices can be regarded as magnetic islands with plasma 
flow. Quasistatic MHD and hydrodynamic cat eyes can be recovered 
as particular cases. Also, it may be noted that for flows non parallel to 
the magnetic field, the electric field is perpendicular to the magnetic 
surfaces [Eq. (9)]. 

3. In fact, (lT2"j) and (fTBl hold for a rather large set of equilibria because 
the functions p(u), $(«), M{u) and one out of B z {u) and P s {u) remain 
free. 

We will further consider a subset of steady sates by assigning the free 
functions P s , B z and M as 

P>) = /?^exp(-2u) + ^, (14) 

B 2 z {u) = (1 - - 6 2 ) exp (-2u) + 5 2 , (15) 
M = M exp(— 2nw) = M (cosh y — e cos x)~ 2n , n > 0. (16) 
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Choice (I16p yields a peaked M 2 -profile along y with |M | being the maxi- 
mum absolute value at x = y = 0. The profile becomes steeper as n takes 
larger positive values, thus increasing the shear of M in relation to the ve- 
locity shear. Henceforth, profiles will refer to the y-axis. The parameter B z q 
represents the external axial magnetic field, 

_ P s (e = y = 0) 
P B 2 z0 /2 

and (3f = P s0 / (B 2 /2), where P s0 = const. Note that ft has been introduced 
in (fT4"|) and (fT5l) in such a way that ffT5j) is automatically satisfied. The other 
parameter ftj in (]14p yields force-free quasistatic equilibria when ft = 0. For 
ft 7^ 0, we set ftf = in order that P s vanishes for y — > oo; thus, only 
one of the parameters ft and /?/ is finite in connection with peaked and flat 
P s -profiles, respectively. For flat P s -profiles, to guarantee positiveness of the 
pressure for ftf > 0, Eq. (j^D is modified to: 

„ ftf M 2 .„ l2 M n 2 

P= — Vn H 5 (17) 

2 2(1 — M 2 ) 2(1 — Mq) ' 1 ' 

The parameters M , and n are free together with L, g , e, P^o, and ft or 
It is recalled that dimensionless quantities are employed and therefore 
B z0 = 1. Also, the reference quantities L and £> not appearing explicitly in 
the equations can arbitrarily be defined as the vortex length (along the x-axis) 
and the density at x — y — 0. Because of the many free parameters, there 
is a variety of steady states including extensions of equilibria employed as 
initial states in reconnection studies (see for example Ref. p2]). An example 
concerns the one-dimensional, force free quasistatic equilibrium recovered for 
ft — e — 0. In the presence of flow and e ^ this equilibrium becomes two- 
dimensional with hollow pressure profile (Eq. (fl7|) ; see also Fig. 4b). Also, in 
this case both current density and velocity have all three components finite. 
Another example for ft = 1 is an equilibrium with B z = B z q, axial current 
density and three-component velocity. For vanishing flow and e = this 
reduces to the Harris sheet equilibrium [T5] . 

We have examined the pressure by Mathematica 6 within broad regions 
of the free parameters, i.e., < e < 1, < ft < 0.9, < M < 0.9 and 
< n < 15. Note that, because of the flow term in (jSJ) the pressure for 
certain parametric values can become negative. Thus, particular care has 
been taken in getting everywhere physically acceptable pressure. For two 
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dimensional equilibria, it turns out that the flow has strong impact on the 
isobaric surfaces by creating "pressure islands" within the cat eyes. This is 
shown in Fig. [31 Also P-profiles are presented in Fig. HI As can be seen 
in Fig. 3 pressure islands appear even for parametric values of experimental 
fusion concern (/3 = 0.02, Mq = 0.02). Since for linear equilibria the flow 
impact on the pressure is weak, it is the nonlinearity here which should play 
an important role. Also, as will be discussed in Sec. Ill, the formation of 
pressure islands may be related to appreciable stabilizing effects of the flow. 

II. Stabilizing effects of the flow 

The stability of the equilibria described by (|T2|) and (|T4l) - (fT6l) is now 
examined by applying a recent sufficient condition [5] . This condition states 
that a general steady state of a plasma of constant density and incompressible 
flow parallel to B is linearly stable to small three-dimensional perturbations 
if the flow is sub-Alfvenic (M 2 < 1) and A > 0, where A is given by Eq. 
(20) of Ref. [5]. Consequently, we restrict the study to parallel flows and 
set g — 1. First it is noted that on the basis of Mercier expansions it turns 
out that the condition is never satisfied in the vicinity of the magnetic axis 
(A < 0) [H]. This holds for generic two-dimensional equilibria irrespective 
of the geometry. Also, for the pressure (fTTl) . the quantity A is independent 
of j3f, as may be expected on physical grounds, because A contains dP s /du 
and not P s itself. In the w-space for translationally symmetric equilibria, A 
assumes the form 



.4 



-S 2 |(jxVu)-(B- V)Vu 
(l - M 
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~ 1/2 Vu- — + g(l-M 2 Y l \Vu\ 2 



M' 



|Vu| 



(18) 



where 



-1/2 



'dP, dM 2 B< 



du 



du 2 



(19) 



and B and j as given by (6) and (7). To calculate A analytically for the 
equilibria under consideration we developed a code in Mathematica 6. The 
expressions obtained for both peaked and flat P^-profiles being lengthy are 
not given explicitly here except for the case of quasistatic equilibria [Eq. (BOl 
below]. The calculations led to the following conclusions. 
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1. For quasistatic equilibria (M = 0) the quantity A assumes the concise 
form 

^ e(l — e 2 ) [ecosh(y) sin(x) 2 + cos(x) sinh(|/) 2 ] 

[e cos(x) — cosh(y)] 5 

Note that A becomes independent of (3 and B z q. The condition is 
nowhere satisfied in the island region except for one dimensional con- 
figurations (e = 0), point vortices (e = 1), the magnetic axes, the 
x-points and for y — > oo for which A = 0. A profile of A is given in 
Fig. B 

2. The flow results in the formation of a stable region close to the magnetic 
axis in the location of pressure islands. An example shown the sign of 
A on the poloidal plane is presented in Fig. 6a. The red colored regions 
are stable (A > 0), while in the blue colored region it holds A < 0. The 
whole area of Fig. 6a becomes blue colored when Mq = 0. 

3. The stable region broadens when the parameters Mq, n and e take larger 
values as can be seen in Figs. 7a, 7b and 7c, respectively. Note the sen- 
sitiveness of A in the region of the stable window to the small variation 
of these parameters possibly related to the nonlinearity; in particular, 
e appears in the argument of the cat-eyes solution (fT2j) . These results 
hold for both peaked- and flat-P s equilibrium profiles. Unlikely, the 
stable region is rather insensitive to the variation of (3. An example is 
given in Fig. 7d, where the stable window persists (just getting slightly 
smaller) when (3 is increased by an order of magnitude (from 0.02 to 
0.2). Also, for point vortices (e = 1) A becomes independent of (3 
irrespective of the value of Mq. 

4. Although for Mq = the vacuum magnetic field B z q has no impact 
on A [Eq. (1201) ]. in combination with the flow, B z0 can enhance the 
stable region. An example of this synergetic effect is shown in Fig. 8a. 
Another example of such a strong synergism can be seen in Fig. 8b 
for a two-dimensional Harris-type equilibrium {(3 = 1, e ^ 0). In this 
case, while the flow itself can not make A positive, together with B z0 
it results in the formation of the stable window. 

V. Summary and Conclusions 
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We have extended the "cat-eyes" solution of the hydrodynamic equilib- 
rium equation to cover MHD magnetically confined plasmas with incompress- 
ible flow. The extension was accomplished smoothly because the pertinent 
generalized Grad-Shafranov equation can be transformed to a form identical 
with that of the hydrodynamic equation [Eq. (|1J)]. Velocity, magnetic field 
and current density of the extended equilibrium share the same surfaces; 
therefore, the vortices can be viewed as magnetic islands with flow with the 
magnetic field not affecting the vortex shape. Also, to be compatible with the 
cat-eyes solution, the axial magnetic field, B z (u), and the quasistatic pres- 
sure, P s (u), must satisfy relation (Tl3l) . The equilibrium is generic enough 
because four surface quantities, i.e. the density, the electrostatic potential, 
the poloidal Alfven Mach function [M(«)] and either B z (u) or P s (u) remain 
free. Generalized Harris or force free-type equilibria can be derived as partic- 
ular cases. Furthermore, the flow caused departure of the pressure surfaces 
from the magnetic surfaces has been examined by assigning the functions 
B z (u), P s (u) and M{u) [Eqs. fll4til6p ]. The equilibrium has the following 
seven free parameters: the island length (L), the density g on the island 
axis, the external axial magnetic field (B z0 ), a parameter e determining the 
island size, a local ratio of the thermal pressure to the magnetic pressure 
((3 or j3f in connection with peaked and flat profiles of P s , respectively), the 
Mach number M on the island axis, and a velocity-shear-related parameter 
n. It turns out that, unlike to linear equilibria, the flow strongly affects the 
pressure surface topology by forming pressure islands on the poloidal plane 
within the cat eyes, even for flows of laboratory fusion concern. 

For parallel flows and constant density, the linear stability of the equilib- 
ria constructed has been examined by means of a recent sufficient condition 
guaranteeing stability when the flow is sub-Alfvenic and an equilibrium de- 
pendent quantity A [Eq. 0181) ] is nonnegative. By symbolic computation of 
A for a broad variation of the parameters e, /3, Mq and n, we came to the 
following conclusions. The flow can result in the formation of a stable region, 
close to the magnetic axis in the location of pressure islands, thus indicating 
a correlation of stabilization with nonlinearity. The stable region can appear 
for fusion relevant values of M on the order of 0.01 when the velocity shear 
becomes appropriately large (n « 10), enhances as n becomes larger and 
persists for a large variation of (3 (from 0.02 to 0.2). Also, the broader the 
stable region the larger is the island size (larger e). A combination of velocity 
and B z0 can have synergetic stabilizing effects by enlarging the stable region. 

In conclusion, the present study has shown significant stabilizing effects 
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of the flow and flow shear in connection with nonlinearity and formation 
of equilibrium pressure islands. The study can be extended to several di- 
rections. Firstly, since four surface functions remain free in equilibrium to- 
gether with many free parameters, there may be a possibility of stability 
optimization. Secondly, the problem could be examined in cylindrical and 
axisymmetric geometries in connection with the magnetic field curvature and 
toroidicity. Note that in the presence of toroidicity non parallel flows have a 
stronger impact on equilibrium because, in addition to the pressure, they re- 
sult in a deviation of the current density surfaces from the magnetic surfaces. 
Although in non-plane geometries nonlinear solutions in general should be 
constructed numerically, it is interesting to pursue analytic translationally 
symmetric solutions in cylindrical geometry as a next step to the cat-eyes 
solution. At last, it is recalled that the search for necessary and sufficient 
stability conditions with flow remains a tough problem as already known in 
the framework of hydrodynamics. 
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Figure captions 



Fig. Q] u-lines of the MHD cat-eyes solution (QJ) for e = 0.2 and (3 = 0.02 
as intersections of the magnetic surfaces with the poloidal plane. 

Fig. [2J Profile of the quasistatic pressure function P s [Eq. (1141) ] along the 
y— axis. 

Fig. [3j Pressure islands in connection with Eqs. ([5]) [Fig. 3a] and (jTTl) 
[Fig. 3b]. The curves represent pressure lines on the poloidal plane. In the 
absence of flow the lines of Fig. 3a coincide with the w-lines of Fig. [T] while 
the equilibrium of Fig. 3b becomes force-free. 

Fig. HJ Pressure profiles along the y-axis respective to the pressure-island 
configurations 3a and 3b. For vanishing flow the profiles 4a and 4b become 
peaked and flat, respectively. 

Fig. Profile of the quantity A [Eq. (FIB!) ] associated with the sufficient 
condition for linear stability for a quasistatic equilibrium (M = 0). Except 
for the marginally stable points y = and y — > oo the condition is nowhere 
else satisfied. 

Fig. O Stabilization effect of flow: In the presence of flow the red colored 
stable regions appear in the diagram 6a where A > 0. The respective stable 
window can be seen in the profile of A in 6b. 

Fig. [7J Impact of the flow (7a), flow shear (7b), cat-eyes size (7c) and 
thermal pressure (7d) in connection with a variation of the parameters Mo, 
n and e, and f3, respectively, on the flow caused stable window associated 
with A > for the equilibrium of Fig. 3a. 

Fig. Combined stabilization effect of flow and B z q\ The curve 8a indicates 
a stabilizing synergism of B z q and flow for the equilibrium of Fig. 3a. A 
stronger synergism of this kind is shown in Fig. 8b pertaining to a two- 
dimensional Harris-type equilibrium. 
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List of Figures 




Figure 1: u-lines of the MHD cat-eyes solution ffl2|) for e = 0.2 and (3 = 0.02 
as intersections of the magnetic surfaces with the poloidal plane. 

e=0.2, p=0.02,p f =0 



0.014 
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0.008 
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0.002 



Figure 2: Profile of the quasistatic pressure function P s [Eq. (TH|) ] along the 
y— axis. 
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3a) e=0.2, y6=0.02, M () =0.02, n=8 3b) e=0.2, ^=0.02, M =0.02, n=8 

y y 




-4 -2 2 4 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



Figure 3: Pressure islands in connection with Eqs. (J5]) [Fig. 3a] and (]17[) 
[Fig. 3b]. The curves represent pressure lines on the poloidal plane. In the 
absence of flow the lines of Fig. 3a coincide with the w-lines of Fig. [1] while 
the equilibrium of Fig. 3b becomes force-free. 



4a) e=0.2,/?=0.02, M =0.02, n=8 4b) e=0.2, ^=0.02, M =0.02, n=8 

P P 




Figure 4: Pressure profiles along the ?/-axis respective to the pressure-island 
configurations 3a and 3b. For vanishing flow the profiles 4a and 4b become 
peaked and flat, respectively. 
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e=0.2, /?=0.02, M =0.02, n=8 

A 




Figure 5: Profile of the quantity A [Eq. ([TBI) ] associated with the sufficient 
condition for linear stability for a quasistatic equilibrium (M = 0). Except 
for the marginally stable points y = and y — ► oo the condition is nowhere 
else satisfied. 



6a) e=0.2, ,3=0.02, M =0.02, n=8, B z0 =0 

y 




Figure 6: Stabilization effect of flow. In the presence of flow the red colored 
stable regions appear in the diagram 6a where A > 0. The respective stable 
window can be seen in the profile of A in 6b. 
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7a) e=0.2,y3=0.02, n=8, 5 z0 =0 



7b) e=0.2, /?=0.02, M =0.02, B z0 =0 



A 



n = 8.5 
n = 8.0 




7c) ,3=0.02, M () =0.02, n=8, S z0 =0 



e = 0.21 
e = 0.20 



7d) e=0.2, M () =0.02, n=8, 5 z0 =0 





Figure 7: Impact of the flow (7a), flow shear (7b), cat-eyes size (7c) and 
thermal pressure (7d) in connection with a variation of the parameters Mo, 
n and e, and f3, respectively, on the flow caused stable window associated 
with A > for the equilibrium of Fig. 3a. 
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8a) e=0.2,y3=0.02, M =0.02, n=8 



8b) e=0.2,yS=l,M () =0.2, n=3 




A 

0.15 r 




Figure 8: Combined stabilization effect of flow and B z0 : The curve 8a in- 
dicates a stabilizing synergism of B zQ and flow for the equilibrium of Fig. 
3a. A stronger synergism of this kind is shown in Fig. 8b pertaining to a 
two-dimensional Harris-type equilibrium. 
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